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By  controlling  the  thermal  transport  of  VO2  nano-gap  junctions  using  device  geometry,  contact  material,  and  applied 
voltage  waveforms,  the  electronically  induced  insulator-metal  phase  transition  is  investigated  in  the  adiabatic  heating 
and  transient  carrier  injection  regimes.  With  a  gradual  ramping  of  an  applied  voltage  on  a  microsecond  time  scale, 
the  transition  electric  field  threshold  can  be  directly  reduced  by  the  Joule  heating.  With  an  abrupt  applied  voltage,  the 
transition  threshold  is  initiated  by  carriers  injected  within  the  first  tens  of  nanoseconds,  but  the  complete  insulator-metal 
phase  transition  is  limited  by  thermal  redistribution  times  to  hundreds  of  nanoseconds. 


The  dynamics  of  the  insulator-metal  phase  transition  (IMT) 
in  the  correlated  electron  material  vanadium  dioxide  (VO2)1’2 
have  attracted  much  attention3,4  due  to  potential  application 
of  VO2  in  transistors5,6  and  opto-electronic7,8  switches.  A  re¬ 
versible  IMT  through  a  structural  reconfiguration  occurs  in 
VO2  when  the  local  temperature  rises  beyond  the  insulator- 
metal  transition  temperature,  Tj^m  —  340  K9.  This  pro¬ 
cess  takes  place  on  the  time  scale  of  10_7-10_5s  for  nano¬ 
scale  devices5.  Alternatively,  the  IMT  due  to  an  electronic 
reconfiguration  takes  place  on  time  scales  of  10_12-10-8  s10 
due  to  a  perturbation  in  the  electronic  structure11,  accumu¬ 
lated  surface  charge12,  external  strain13,  applied  electric  fields 
or  voltages14-18,  and  currents19.  An  electronically  induced 
IMT  produced  by  a  current  or  applied  voltages  is  typically 
accompanied  by  a  local  temperature  increase  in  the  VO2  from 
Joule  heating20,21.  The  temperature  rise  modifies  the  tempo¬ 
ral  dynamics  of  the  phase  transition  by  inducing  a  secondary 
structural  phase  transition21  when  the  temperature  exceeds 
Although  thermal  effects  can  be  reduced  by  limit¬ 
ing  the  current22  or  increasing  the  heat  dissipation  rate  of  the 
devices23,  the  role  of  Joule  heating  during  an  electronically 
induced  IMT  remains  unclear5,24. 

In  this  Letter,  we  report  the  effects  of  heating  on  the  elec¬ 
tronically  induced  IMT  in  VO2  nano-gap  junctions.  The  ther¬ 
mal  transport  was  controlled  by  the  applied  voltage  waveform, 
device  geometry,  and  metal  contact.  We  studied  the  IMT  in 
the  adiabatic  and  transient  regimes.  The  dependence  of  the 
IMT  threshold  voltage  on  the  duration  of  the  applied  voltage 
pulses  and  device  parameters  can  be  explained  by  thermal  ac¬ 
cumulation  and  dissipation  rates.  The  IMT  could  not  occur  in 
the  absence  of  injected  carriers,  so  carrier  effects  are  critical 
in  the  electronically  induced  phase  transition  of  VO2. 

The  schematic  and  scanning  electron  micrograph  (SEM)  of 
the  typical  device  are  presented  in  Figs.  1(a)- 1(b)  respectively. 
A  VO2  film  with  a  thickness  of  h  =  100  nm  was  deposited  us¬ 
ing  radio-frequency  magnetron  sputtering  of  a  vanadium  tar¬ 
get  on  2  pm  thermally  grown  layer  of  silica  (Si02)  on  a  silicon 
(Si)  substrate.  The  VO2  junctions  were  formed  by  a  combi¬ 
nation  of  electron  beam  (e-beam)  lithography,  metal  evapora¬ 
tion,  and  lift-off  to  produce  two  metal  contacts  separated  by  a 


FIG.  1.  (a)  The  top  and  side  view  schematic  diagram  of  a  VO2  lat¬ 
eral  junction,  (b)  SEM  of  the  device  and  the  magnified  image  of  the 
gap.  (c)  The  simplified  model  for  the  heat  dissipation  and  generation 
mechanisms  in  the  device,  (d)  The  simplified  circuit  model  of  the 
device  and  the  experiment. 


gap  of  length  L.  These  contacts  had  a  width  of  W  =  10  pm 
and  a  thickness  of  t  =  100  nm.  The  contacts  were  either  gold 
(Au)  or  palladium  (Pd). 

The  phase  transition  dynamics  are  limited  by  the  RC  time 
constant  of  the  circuit,  trc,  the  thermal  extraction  time  of 
the  junction,  77,  and  the  forward  electronic  phase  transitions 
time  constant,  77.  An  independent  measurement  of  te  and  tt 
of  the  junctions  is  difficult,  since  the  electronic  phase  transi¬ 
tion  is  typically  accompanied  with  a  temperature  increase.  As 
the  temperature  increases,  the  average  electric  field  at  the  on¬ 
set  of  the  phase  transition,  decreases25  according  to 

Ej^m  oc  (1  —  (T/Tc))2.  However,  by  independently  vary¬ 
ing  tt,  we  can  evaluate  the  roles  of  the  thermal  and  electronic 
effects,  as  well  as  their  interaction  in  the  phase  transition. 

The  heat  transfer  channels  in  the  VO2  junction  are  illus¬ 
trated  in  Fig.  1(c).  The  total  heat  dissipation  rate  of  the  junc¬ 
tion  is  the  sum  of  the  effective  heat  extraction  rates  of  the  sur¬ 
rounding  media,  7^,  where  x  is  either  air,  VO2,  SiC>2  or  the 
metal  contact.  Since  the  thermal  conductivity  of  the  contact 
metals  is  much  larger  than  that  of  the  VO2,  air,  or  SiC>2,  the 
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FIG.  2.  (a)  Time  profiles  of  a  VbaP  and  Vapp,  when  a  triangular  pulse 
with  £p  =250  qs  was  applied,  (b)  The  measured  as  a  function 

of  L  for  pulses  with  different  tv.  The  measured  (c)  Ei^m  and  (d) 
Vo  as  a  function  of  tp  for  devices  with  Pd  or  Au  contacts. 


total  heat  dissipation  rate  of  the  junction  can  be  approximated 
as  7Metai-  The  heat  generation  rate  in  the  junction  is  the  total 
dissipated  electrical  power,  P(t).  As  a  result,  the  simplified 
heat  transfer  equation  for  the  VO2  junctions  can  be  written  as: 

PdCpE=kV^T-^iT+E^,  (1) 

dt  V  eff 

where  is  the  effective  volume  of  the  junction,  and  pv,Cp, 
and  hi  are  the  density,  the  specific  heat  capacity  and  the  ther¬ 
mal  conductivity  of  VO2  respectively.  This  equation  shows 
that  7Metai  and  P(t )  changes  the  thermal  transport  time  of  the 
system  and  allows  for  the  controlled  study  of  77  and  te.  In 
the  experiments  to  follow,  the  contacts  were  either  Au  or  Pd, 
and  we  controlled  P(t)  by  the  applied  voltage  waveform.  The 
thermal  conductivities  of  Au  and  Pd  are  about  315  W/(m-K) 
and  70  W/(m-K),  respectively. 

To  measure  the  time  constants  of  the  VO2  devices,  an  arbi¬ 
trary  function  generator  was  used  to  apply  a  voltage,  Vapp, 
across  the  device  contacts  and  a  series  resistor  RL  =  550 
fi  .  The  simplified  circuit  model  of  the  device  and  the  experi¬ 
ment  is  shown  in  Fig.  1(d).  The  voltage  drop  across  the  gap, 
Voap,  and  Fapp  were  monitored  using  an  oscilloscope.  The 
measured  trc  for  the  setup  was  «  40  =b  10  ns  and  remained 
constant  within  the  experimental  error  range  for  devices  with 
different  values  of  gap  lengths  and  contact  metals. 

We  investigated  the  adiabatic  Joule  heating  regime  by  ap¬ 
plying  a  periodic  train  of  triangular  voltage  pulses,  which 
gradually  increased  I  and  Vbap  and  allowed  sufficient  time  for 
Joule  heating  to  occur.  The  rise  time  of  the  pulses,  £p,  was 
at  least  5  qs.  The  period  was  2  s  to  prevent  cumulative  heat¬ 
ing  from  successive  pulses.  Figure  2(a)  shows  an  example 
of  the  measured  Vbap  and  Vapp.  The  phase  transition  mani¬ 
fested  in  a  drop  in  Vbap  at  the  forward  phase  transition  voltage 
Vi^m,  and  an  increase  in  Voap  at  the  reverse  transition  volt- 
age,  VM^i26. 


To  extract  the  electric  field  for  the  phase  transition,  we  mea¬ 
sured  Vi^m  for  devices  with  different  L  and  took  the  slopes 
of  to  be  the  “average”  transition  electric  field,  Ej^m-  The 
voltage  offset,  Vo,  taken  to  be  the  value  of  Vi^m  as  L  -A  0, 
is  due  to  the  finite  contact  resistance  and  work  function  mis¬ 
match  between  VO2  and  the  contact  material.  As  an  example, 
Figure  2(b)  shows  the  measured  Vj^m  as  a  function  of  L  for 
devices  with  Pd  contacts  and  for  tv  between  5  qs  and  1  s.  For 
all  device  lengths,  shorter  voltage  pulses  required  higher  val¬ 
ues  of  Vi->m,  which  meant  more  energy  had  to  be  deposited 
to  the  device  to  initiate  the  transition.  This  could  be  caused 
by  a  reduction  of  thermal  accumulation  (which  would  have 
decreased  V/_^m25)  and/or  a  limitation  due  to  te. 

To  investigate  thermal  influence  on  Ej^m,  we  extracted 
Ei^m  as  a  function  of  tv  for  devices  with  Au  or  Pd  metal 
contacts  in  Fig.  2(c).  The  solid  fines  are  fits  to  guide  the  eye. 
For  both  types  of  samples,  Ei^m  was  constant  when  tp  >  1 
ms.  As  tp  decreased,  Ej^m  increased  and  saturated  at  around 
55  kV/cm  for  the  Au  samples  when  tp  <  10-5s.  Ei^m  did 
not  saturate  for  the  Pd  samples  for  the  measured  values  of  tp. 
We  define  a  characteristic  pulse  width,  tT,  as  the  value  of  tp 
that  resulted  in  a  50%  increase  in  Ei^m  from  its  value  at 
tp  =  Is.  As  shown  in  Fig.  2(c),  tju  ^3x  10-5s  for  the  Au 
samples  and  tjd  «  5  x  10-6s  for  the  Pd  samples. 

The  dependence  of  Ej^m  on  tp  and  the  contact  metal  can 
be  explained  by  thermal  transport.  When  tp  tT,  Joule  heat¬ 
ing  raised  the  temperature  in  the  VO2  junction,  which  reduced 
Ei^m-  As  tp  decreased,  heating  had  a  lesser  contribution  and 
consequently  Ei^m  increased.  Since  Au  has  a  higher  thermal 
conductivity  than  Pd  (which  would  reduce  the  temperature  of 
the  VO2  in  the  Au  samples),  Ei^m  for  Au  was  larger  than 
in  the  Pd  samples  at  all  values  of  tp.  Ej^m  saturated  when 
tv  tj,  which  was  approximately  the  value  of  Ej^m  when 
no  thermal  effects  were  present.  In  the  Au  samples,  the  satu¬ 
ration  of  E^m  to  about  55  kV/cm  when  tp  <C  tr  occurred 
when  almost  no  thermal  effects  were  present.  This  field  value 
agrees  well  with  measurements  of  vertical  VO2  junctions,  in 
which  heat  accumulation  is  minimized  owing  to  large  device 
cross-sections27. 

The  temperature  change  in  the  VO2  due  to  the  applied  volt¬ 
age  further  affects  the  carrier  density  in  the  devices.  Fig.  2(d) 
shows  the  offset  voltage,  Vo,  vs.  tp.  Both  the  carrier  density28 
and  the  work  function29  of  VO2  increase  across  the  phase 
transition  with  a  similar  temperature  dependence  as  Ej^m25- 
However,  while  a  higher  carrier  density  reduces  Vo,  a  higher 
work  function  increases  Vo.  The  carrier  density  effects  tend 
to  dominate  over  the  work  function  changes,  because  the  car¬ 
rier  density  can  change  by  about  3  orders  of  magnitude  across 
the  IMT,  while  the  work  function  changes  by  less  than  3%29. 
The  measured  resistivity  change  of  the  junctions  was  slightly 
larger  than  2  orders  of  magnitude23 .  Therefore,  the  increase 
in  Vo  when  tp  <  tT,  in  a  similar  way  as  Ej^m  in  Fig.  2(c), 
suggests  the  carrier  density  decreased  as  the  voltage  pulses 
became  shorter  and  the  role  of  heating  was  reduced. 

To  study  the  IMT  without  thermal  effects,  we  applied 
square- wave  pulses  with  rise  time  tp  <C  tj,  which  was  lim¬ 
ited  to  40  ns  due  to  the  experimental  setup.  The  pulses  had  a 
duration  of  1  qs  and  a  period  of  2  s.  The  amplitude  of  Vapp, 
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FIG.  3.  (a)  The  measured  Vapp  and  VbaP  when  a  square- wave  pulse 
with  a  duration  of  1  ji s  was  applied  across  the  devices,  (b)  The  mea¬ 
sured  resistance  of  the  device  as  VbaP  increased.  The  measured  (c)  n 
and  (d)  T2  of  the  devices  vs.  the  applied  voltage  amplitude. 


^App,  was  varied  and  the  time  profile  of  Vbap  was  measured. 
Figure  3(a)  shows  an  example  of  the  measured  Vapp  and  Vbap, 
and  i?App  is  indicated.  As  before,  the  drop  in  Vbap  indicated  the 
start  of  the  phase  transition.  The  value  of  Ei^m  was  deter¬ 
mined  using  the  smallest  value  of  vavv  that  triggered  the  IMT 
for  each  L  and  was  52  db  5  kV/cm  for  Pd  and  55  db  8  kV/cm 
for  Au  devices.  These  values  agreed  well  with  the  saturation 
values  of  Ej^m  in  Fig.  2(c). 

Figure  3(b)  shows  the  measured  values  of  i?Gap  as  t>app  in¬ 
creased  for  a  device  with  L  =  800  nm  and  Pd  contacts.  i?Gap 
was  constant  when  Vbap  <  Vi^m-  Increasing  t;App  such  that 
Vbap  >  Vi^m  resulted  in  a  rapid  drop  in  i?Gap  at  a  time  of 
t i  after  the  rise  of  the  pulse.  Rq ap  continued  to  decrease  and 
dropped  below  10  dB  of  its  initial  value  at  r 2  after  the  rise 
of  the  pulse.  Both  t\  and  r 2  decreased  with  increasing  vApp. 
Figures  3(c)- 3(d)  show  the  measured  t\  and  72  as  a  function 
of  ^APP  for  devices  with  two  different  values  of  L  and  contact 
metals.  As  vavv  increased,  t\  for  the  L  —  500  nm  Pd  and  Au 
devices  and  L  =  800  nm  Pd  devices  decreased  and  saturated 
at  approximately  trc.  This  suggests  t\  is  due  to  the  combined 
effects  of  trc  and  tr  with  tr  <  trc-  t\  did  not  saturate  for 
the  L  =  800  nm  Au  device  for  the  voltages  tested.  r2  also 
decreased  with  increasing  7>App  at  a  lower  rate  than  ti,  con¬ 
sistent  with  T2  being  caused  by  Joule  heating  of  the  VO2.  As 
expected,  longer  channels  required  more  time  for  the  heat  to 
redistribute.  For  the  same  device  length  and  the  same  voltage 
difference  from  the  phase  transition  threshold,  r2  for  Au  was 
higher  than  that  for  Pd  by  up  to  a  factor  of  about  2,  because 
the  local  temperature  of  the  Au  samples  would  be  lower  owing 
to  the  higher  thermal  conductivity  of  Au.  As  t^pp,  or  equiv¬ 
alently,  the  applied  electrical  power,  increased,  r 2  decreased 
because  of  the  faster  rate  of  heating  as  described  by  Eq.  1 . 

These  measurements  confirm  that  the  phase  transition  is  ini¬ 
tiated  electronically  and  occurs  on  a  timescale  limited  by  the 
RC  time  constant  of  the  measurement  setup,  and  is  followed 
by  a  thermal  transition  that  is  limited  by  the  thermal  dissipa¬ 


FIG.  4.  (a)  The  top  and  side  view  schematics  of  a  VO2  lateral  junc¬ 
tion  with  a  dielectric  insulating  layer,  (b)  SEM  image  of  the  device  at 
the  gap.  The  contrast  of  the  ZrC>2  has  been  enhanced  using  an  X-ray 
detector,  (c)  The  TV  characteristics  of  device  with  L  —  500  nm.  The 
measured  data  are  marked  by  the  dots,  and  the  dashed  lines  connect 
the  data  across  the  phase  transition,  (d)  The  measured  Vi^m  vs.  L 
for  devices  with  and  without  the  ZrC>2  layer. 


tion  of  the  junction.  To  suppress  the  thermal  effects,  it  is  es¬ 
sential  to  reduce  I  by  electrically  insulating  the  VO2  from  the 
metallic  contacts.  One  approach  is  to  use  a  thin  high-k  dielec¬ 
tric  between  the  VO2  and  the  metallic  contacts.  Figure  4(a) 
shows  the  schematic  of  the  device.  A  7  nm  thick  layer  of  zir¬ 
conium  dioxide  (Zr02)  was  deposited  using  atomic  layer  de¬ 
position  and  patterned  using  a  combination  of  electron  beam 
lithography  and  lift-off.  The  metal  contacts  were  formed  as 
before  using  Pd  and  had  the  width  of  TV  =  10  gm  and  thick¬ 
ness  of  100  nm.  Figure  4(b)  is  an  X-ray  enhanced  SEM  image 
of  the  device  near  the  gap  showing  that  the  thin  dielectric  layer 
that  covers  the  VO2  surface. 

To  observe  the  phase  transition,  we  measured  Vbap  and  the 
current,  /,  that  leaked  through  the  ZrC>2  and  VO2  layers.  Fig¬ 
ure  4(c)  shows  the  measured  IV  characteristics  of  two  junc¬ 
tions  with  (green)  and  without  (black)  the  ZrC>2  layer.  Both 
devices  had  the  same  L  =  500  nm.  The  current,  /,  initially  in¬ 
creased  with  increasing  Vbap  and  the  phase  transition  resulted 
in  a  drop  in  Vbap-  As  expected,  the  ZrC>2  layer  reduced  I  by 
an  order  of  magnitude  at  any  value  of  Vbap-  However,  the 
ZrC>2  layer  also  increased  Vi^m  by  a  factor  of  2.5,  which 
resulted  in  a  value  of  I  that  was  similar  to  devices  without 
ZrC>2.  Figure  4(d)  shows  the  measured  Vi^m  vs.  L  for  de¬ 
vices  with  (green)  and  without  (black)  the  ZrC>2  layer.  Vi^m 
did  not  change  with  L  for  devices  with  the  ZrC>2  layer,  indicat¬ 
ing  that  the  phase  transition  occurred  due  to  the  electric  field 
breakdown  of  the  ZrC>2  that  allowed  sufficient  current  to  pass 
through  the  device  before  the  phase  transition  occurred. 

As  another  approach  to  limit  the  current,  we  fabricated  VO2 
nano-  and  micro-wires  placed  between  a  capacitor  with  an  air 
gap.  The  VO2  wires  were  fabricated  using  e-beam  lithogra¬ 
phy  and  subsequent  reactive  ion  etching.  The  metal  contacts 
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FIG.  5.  (a)  SEM  images  of  the  device,  VO2  nanowire  in  the  capacitor 
gap,  and  metal  contact,  (b)  The  normalized  resistivity  of  the  wire,  (c) 
The  wire  resistivity  as  a  function  of  Vbap  at  different  temperatures. 


were  defined  using  e-beam  lithography,  thermal  evaporation 
of  400  nm  of  Au,  and  a  subsequent  lift-off  step.  Figure  5(a) 
shows  the  SEM  images  of  the  fabricated  device  and  the  mag¬ 
nified  images,  showing  the  VO2  nanowire  at  the  gap  and  on 
the  metal  contacts.  Four  metal  pads  contacted  the  VO2  wire 
to  accurately  measure  the  wire  resistivity,  pwire-  A  voltage  of 
Vbap  was  applied  across  the  capacitor,  which  established  an 
electric  field  in  a  section  of  the  VO2  between  the  electrodes. 

Figure  5(b)  shows  the  normalized  resistivity  of  the  wire 
Pwire(^)  / Pwhe(T  =  295  K)  as  a  function  of  T  for  a  wire  with 
a  width  of  500  nm  and  a  thickness  of  100  nm.  When  T  >  335 
K,  the  VO2  underwent  a  phase  transition  and  the  resistivity  of 
the  wires  dropped.  Figure  5(c)  shows  the  normalized  resistiv¬ 
ity  as  Vbap  increased.  The  resistivity  did  not  change  with  Vbap 
up  to  the  breakdown  voltage  of  the  gap,  which  was  about  20  V. 
Below  the  breakdown  voltage,  the  leakage  current  remained 
less  than  10  pA.  The  measurements  were  repeated  at  differ¬ 
ent  temperatures  to  decrease  the  required  Ei^m-  Even  when 
T  «  the  normalized  resistivity  was  unaffected  by 

Vbap-  As  previously  mentioned,  the  required  electric  field  for 
the  phase  transition  of  VO2  scales  as  1/T  and  is  <  5  kV/cm 
close  to  According  to  finite  element  simulations,  the 

average  electric  field  in  these  devices  was  >10  kV/cm.  How¬ 
ever,  these  electric  field  values  did  not  result  in  any  detectable 
changes  in  pWjre  even  at  elevated  temperatures.  If  an  external 
electric  field  could  induce  a  phase  transition  in  these  films,  its 
threshold  would  have  been  >>10  kV/cm,  a  value  which  was 
not  achieved  in  this  experiment.  Consequently,  we  conclude 
that  the  most  plausible  cause  of  the  phase  transition  in  these 
films  was  carrier  injection. 

In  summary,  we  have  shown  the  influence  of  Joule  heat¬ 
ing  on  the  threshold  electric  field,  carrier  densities,  and  phase 
transition  time  scales  of  the  IMT  in  micro-  and  nano- scale  two 
terminal  VO2  junctions.  The  phase  transition  was  initiated  by 
carrier  injection  (not  electric  field).  The  measured  thermal  dis¬ 


sipation  time  was  10“ 7  s,  while  the  measured  electronic  tran¬ 
sition  time  was  limited  by  the  RC  time  constant  of  the  setup 
to  40  ns.  The  thermal  effects  can  be  suppressed  using  high 
frequency  circuits  for  current  regulation  after  the  phase  tran¬ 
sition.  These  circuits  are  typically  used  for  current  quenching 
in  avalanche  photodiodes30  and  can  be  potentially  adapted  to 
VO2  devices. 
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